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Abstract 

This paper shows that the Camassa-Holm (CH) spectral problem yields two 
different integrable hierarchies of nonlinear evolution equations (NLEEs), one 
is of negative order CH hierachy while the other one is of positive order CH 
hierarchy. The two CH hierarchies possess the zero curvature representations 
through solving a key matrix equation. We see that the well-known CH equation 
is included in the negative order CH hierarchy while the Dym type equation is 
included in the positive order CH hierarchy. Furthermore, under two constraint 
conditions between the potentials and the eigenfunctions, the CH spectral prob- 
lem is cast in: 

1. a new Neumann- like iV-dimensional system when it is restricted into a 
symplectic submanifold of M. 2N which is proven to be integrable by using 
the Dirac-Poisson bracket and the r-matrix process; and 



2. a new Bargmann-like iV-dimensional system when it is considered in the 
whole R 2N which is proven to be integrable by using the standard Poisson 
bracket and the r-matrix process. 

In the paper, we present two 4x4 instead of N x N r-matrix structures. One is 
for the Neumann-like system (not the peaked CH system) related to the negative 
order CH hierarchy, while the other one is for the Bargmann-like system (not the 
peaked CH system, either) related to the positive order hierarchy. The whole 
CH hierarchy (a integro-differential hierarchy, both positive and negative order) 
is shown to have the parametric solutions which obey the corresponding con- 
straint relation. In particular, the CH equation, constrained to some symplectic 
submanifold in R 2N , and the Dym type equation have the parametric solutions. 
Moreover, we see that the kind of parametric solution of the CH equation is 
not gauge equivalent to the peakons. Solving the parametric representation of 
solution on the symplectic submanifold gives a class of new algebro-geometric 
solution of the CH equation. 
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1 Introduction 



The shallow water equation derived by Camassa-Holm (CH) in 1993 [0] is a new in- 
tegrable system. This equation possesses the bi-Hamiltonian structure, Lax pair and 
peakon solutions, and retains higher order terms of derivatives in a small amplitude 
expansion of incompressible Euler's equations for unidirectional motion of waves at the 
free surface under the influence of gravity. In 1995 Calogero || extended the class of 
mechanical system of this type. Later, Ragnisco and Bruschi [^(J and Suris showed 
that the CH equation yields the dynamics of the peakons in terms of an iV-dimensional 
completely integrable Hamiltonian system. Such kind of dynamical system has Lax 



pair and aniVxiV r- matrix structure [20 



Recently, the algebro-geometric solution on the CH equation and the CH hierarchy 
arose much more attraction. This kind of solution for most classical integrable PDEs 
can be obtained by using the inverse spectral transform theory, see Dubrovin 1981 fill , 



Ablowitz and Segur 1981 @, Novikov et al 1984 |7j, Newell 1985 |6|. This is done 



usually by adopting the sprectral technique associated with the corresponding PDE. 
Alber and Fedorov @, ^| studied the stationary and the time-dependent quasi-periodic 
solution for the CH equation and Dym type equation through using the method of trace 
formula || and Abel mapping and functional analysis on the Riemann surfaces. Later, 
Alber, Camassa, Fedorov, Holm and Marsden || considered the trace formula under the 
nonstandard Abel-Jacobi equations and by introducing new parameters presented the 
so-called weak finite-gap piecewise-smooth solutions of the integrable CH equation and 
Dym type equations. Very recently, Gesztesy and Holden |L4| discussed the algebro- 



geometric solutions for the CH hierarchy using the polynomial recursion formalism and 
the trace formula, and connected a Riccati equation to the Lax pair of the CH equation. 

The present paper is providing another approach to algebro-geometric solutions of 
the CH equation which is constrained to some symplectic submanifold. Our approach 
differs from the ones pursued in Refs. |2|, [| £|, [|, Q and we will outline the differences 
next. Based on the nonlinearization technique |J, we constrain the CH hierarchy 
to some symplectic submanifold and use the constraint between the potentials and 
the eigenf unctions first to give the parametric solution and then to give the algebro- 
geometric solution of the CH equation on the symplectic submanifold. 

The main results of this paper are twofolds. 

• First, we extend the CH equation to the negative order CH hierarchy, which is 
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a hierarchy of integrable integro-differetial equations, through constructing the 
inverse recursion operator. This hierarchy is proven to have Lax pair through 
solving a key matrix equation. The CH spectral problem associated with this 
hierarchy is constrained to a symplectic submanifold and naturally gives a con- 
straint between the spectral function and the potential. Under this constraint, the 
CH spectral problem (linear problem) is nonlinearized as a new iV-dimensional 
canonical Hamiltonian system of Neumann type. This iV-dimensional Neumann- 
like system is not the peaked dynamical system of the CH equation because 
the peakons does not come from the CH spectral problem. The Neumann-like 
CH system is shown integrable by using the so-called Dirac-Poisson brackets on 
the symplectic submanifold in M. 2N and r-matrix process. Here we present a 
4x4 r-matrix structure for the Neumann-like system, which is available to get 
the algebro-geometric solution of the CH equation on this symplectic submani- 
fold. The negative order CH hierarchy is proven to have the parametric solution 
through employing the Neumann-like constraint relation. This parametric so- 
lution does not contain the peakons @, and vice versa. Furthermore, solving 
the parametric representation of solution on the symplectic submanifold gives an 
algebro-geometric solution for the CH equation. We point out that our algebro- 
geometric solution (see Eq. (|3.104j) and Remarks 3 and 4) is different from the 



ones in Refs. 0, 0, f|, |5j, |14|| , and more simpler in the form. 

Second, based on the negative case, we naturally give the positive order CH hier- 
archy by considering the recursion operator. This hierarchy is shown integrable 
also by solving the same key matrix equation. The CH spectral problem, related 
to this positive order CH hierarchy, yields a new integrable iV-dimensional sys- 
tem of Bargmann type (instead of Neumann type) by using the standard Poisson 
bracket and r-matrix procedure in M? N . A 4 x 4 r-matrix structure is also pre- 
sented for the Bargmann- like system (not peaked CH system, either), which is 
available to get the parametric solution of a Dym type equation contained in 
the positive order CH hierarchy. This hierarchy also possesses the parametric 
solution through using the Bargmann constraint relation. 



Roughly speaking, our method works in the following steps (also see 

- Start from the spectral problem. 

- Find some constraint condition between the potentials and the eigenfunctions. 
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Here, for the negative CH hierarchy, we restrict it to a symplectic manifold in ~R 2N , but 
for the positive CH hierarchy, we will have the constraint condition in the whole M. 2N . 

- Prove the constrained SPECTRAL PROBLEM is fmite-dimenaional integrable. 
Usually we use Lax matrix and r-matrix procedure. 

- Verify the above constrained potential(s) is (are) a parametric solution of the 
hierarchy. 

- Solve the parametric representation of solution in an explicit form, then give 
the algebro-geometric solutions of the equations on the symplectic manifold. In this 
process, we separate the variables of the Jacobi-Hamiltonian system plU , then construct 
the actional variables and angle-coordinates on the symplectic submanifold, and the 
residues at two infinity points for some composed Riemann-Theta functions give the 
algebro-geometric solutions. 

The whole paper is organized as follows. Next section gives a general structure of the 
zero curvature representations of the all vector fields for a given isospectral problem. 
The key point is to construct a key matrix equation. In section 3, we present the 
negative order CH hierarchy based on the inverse recursion operator. The well-known 
CH equation is included in the negative order hierarchy, and the CH spectral problem 
yields a new Neumann-like system which is constrained to a symplectic manifold. This 
system has the canonical form and is integrable by using the Dirac-Poisson bracket and 
the r-matrix process. Here we obtain a 4 x 4 r-matrix structure for the Neumann-like 
CH system. Furthermore, the whole negative order CH hierarchy constrained on the 
symplectic submanifold has the parametric solution. In particular, the CH equation 
has the parametric solution on this submanifold. Finally we give an algebro-geometric 
solution of the CH equation on this submanifold. In section 4, we deal with the positive 
order integrable CH hierarchy and give a new Bargmann-like integrable system. By the 
use of the similar process to section 3, the CH spectral problem is nonlinearized to be 
an integrable system under a Bargmann constraint. This integrable Bargmann system 
also has an r-matrix structure of 4 x 4. Moreover, the positive order CH hierarchy 
is also shown to have the parametric solution which obeys the Bargmann constraint 
relation. In particular, the Dym type equation in the positive order CH hierarchy has 
the parametric solution. 
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Let us now give some symbols and convention in this paper as follows: 



fg x k f — fkx j k — 0,1,2,..., 
/ ... / fdx, k 1.-2 




ft = §£, fkxt = §3 (k = 0,1,2,...), d = J-, d 1 is the inverse of d, i.e. dd 1 = 



In the following the function m stands for potential, A is assumed to be a spectral 
parameter, and the domain of the spatial variable x is Q which becomes equal to 
(— oo, +oo) or (0, T), while the domain of the time variable is the positive time 
axis IR + = {tfc| t k G R, £fc > 0, k = 0, ±1, ±2, ...}. In the case Q = (—oo, +oo), the 
decaying condition at infinity and in the case Q = (0, T), the periodicity condition for 
the potential function is imposed. 

(M 2Ar ,d]9 A dq) stands for the standard symplectic structure in Euclid space M? N = 
{{p,q)\ V = (Pi, • • • ,Pn), Q = (<?i, ■ ■ ■ , Qn)}, Pj, Qj (j = l,...,N) are pairs of canoni- 
cal coordinates, (•, •) is the standard inner product in ¥L N ; in (M 27V , dpAdq), the Poisson 
bracket of two Hamiltonian functions F, H is defined by 



Ai, Atv are assumed to be N distinct spectral parameters, A = diag(X\, Ajv), and 
^2x2 = diag(l,l). Denote all infinitely times differentiable functions on real field ffi. 
and all integers by C foo (R) and by Z, respectively. 

2 The Camassa-Holm (CH) spectral problem and 
zero curvature representation 

Let us consider the Camassa-Holm (CH) spectral problem 0: 





(1.1) 



i^xx = jip - -mXifj 



(2.1) 
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with the potential function m. 

Eq. ( |2.1|) is apparently equivalent to 





! 1, 



1 



j — ~m\ 



y x = Uy, U = U(m,X) 
where y — (y\, y<i) T = (tp, if) x ) T . Easy to see Eq. fl2.2|) 's spectral gradient 

Xyl 



(2.2) 



$X . , 
VA = — - v ■- 
dm 



(2.3) 



satisfies the following Lenard eigenvalue problem 

K ■ VA = A J ■ VA 
with the pair of Lenard's operators 

K = —d 3 + <9, J = dm + md. 
They yield the recursion operator 

£ = J~ l K = (dm + md)- 1 (d - d 3 ). 

which also has the product form £ = ^m^^d^ 1 m^^(d — d 3 ). 

Apparently, the Gateaux derivative matrix £/*(£) of the spectral matrix U in the 
direction £ e C°°(M) at point m is 



(2.4) 



(2.5) 



(2.6) 



A 



de 



U(m + eO 



e=0 







(2.7) 



which is obviously an injective homomorphism. 

For any given C°°-function G, we construct the following matrix equation with 
respect to V = V{G) 



V x - [U, V] = U*(K -G-XJ-G). 



(2- 



Theorem 1 For the CH spectral problem (\2.3j) and an arbitrary C°° -function G, the 
matrix equation (\2.8[ ) has the following solution 

V = V(G) = x( x 



\G X 



-G 



2&xx — \G + \mXG \G X 



(2.9) 
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Proof: Directly substituting Eqs. ( ET9|) , and ( |2~7D into Eq. ( ^8|) , we can 

complete the proof of this theorem. 



Theorem 2 Let G G #er J = {G G C°°(M) | JG = 0} and G_i G Ker if = {Ge 

C°°(R) | KG = 0}. VFe define the Lenard's sequences as follows 

Gj = C? ■ G = C j+1 ■ G_i, j G Z. (2.10) 

1. the all vector fields X k = J ■ G k , k G Z satisfy the following commutator repre- 
sentation 

V k , x - [U, V k ] = U*(X k ), \/k G Z; (2.11) 

2. the following hierarchy of nonlinear evolution equations 

m tk = X k = J ■ G k , V&GZ, (2.12) 
possesses the zero curvature representation 

U tk -V k , x + [U,V k ] = 0, VkeZ, (2.13) 

where 

and = V(Gj-) gwen 6?/ Eg. fl<?.^ u>z£/i G = Gj. 
Proof: 

1. For = 0, it is obvious. For k < 0, we have 

-l 



v k ,x- [u,v k ] = -Y,(\'j..r- r.Vj)^ i 1 
j=fc 



= U* (j2 K ■ Gj-i\ k ~ j - K ■ GjX^- 1 
\j=k 

= U,{K-G k ^-K-G^\ k ) 
= U^K-Gh-x) 
= U m (X k ). 

For the case of k > 0, it is similar to prove. 
2. Noticing U tk = U*(m tk ), we obtain 

U th - V k , x + [U, V k ] = U*(m tk - X k ). 
The injectiveness of U* implies item 2 holds. 

3 Negative order CH hierarchy, integrable Neumann- 
like system and algebro-geometric solution 

3.1 Negative order CH hierarchy 

Let us first give the negative order hierarchy of the CH spectral problem ( |2.2j ) through 
considering the kernel element of Lenard's operator K. The kernel of operator K has 
the following three seed functions: 

G-i = 1, (3.1) 
G\ = e\ (3.2) 
G 3 _i = e~ x , (3.3) 

whose all possible linear combinations form the whole kernel of K. Let G-\ G Ker K, 
then 

3 

i=i 

where ai = ai(t n ), I = 1,2,3, are three arbitrarily given C°°-functions with respect 
to the time variables t n (n < 0, n e Z), but independent of the spacial variable x. 
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Therefore, G_i directly generates an isospectral (Xt k = 0, k < 0, k e Z) hierarchy of 



(3.5) 



nonlinear evolution equations for the CH spectral problem (2.2) 

m = JC k+1 ■ fc < 0, fc G Z, 



which is called the negative order CH hierarchy because of k < 0. In Eq. ( |3.5| 
the operator J is defined by Eq. ( [2.5D and Cr x is given by 



CT X = K- l J = d- 1 e x d- 1 e- 2x d~ 1 e x {dm + md). 



Here, 



R- 1 = d- 1 e x d- 1 e- 2x d- 1 e x . 
With setting m = u — u xx , we obtain another form of C~ l : 

CT X = u + e x d- 1 e- 2x d~ 1 e x (ud + 2u x + d- 1 m)d. 



(3.6) 



(3.7) 



(3.f 



By Theorem ^, the negative CH hierarchy Q3.5p has the zero curvature representa- 



tion 



^-^ + [[7,^4] 



o, fc < o, k e z, 



(3.9) 
(3.10) 



i.e. 



2/^ 







\-~m\ 



(3.11) 



A; = -1,-2,..., 



where Gj = ■ G_i, j < 0, j 6 Z. Thus, all nonlinear equations in the negative 
CH hierarchy (|3.5|) are integrable. 



Let us now give some special reductions of Eq. ( |3.5|) . 



In the case of a\ = — l,a% = = 0, i.e. G-i = —G^ ± = —1, because G_ 2 
■ G-i = —u, the second equation of Eq. (|3.5| ) reads 
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m t _ 2 = —(dm + md) ■ u, 
i.e. (here noticing m = u — u xx ) 

Ut_ 2 u xx t_„ 4~ 3uu x = 2u x ii xx + uu x 



(3.12) 



(3.13) 



which is exactly the Camassa-Holm equation f7|. According to Eq. ( |3.11|) , the CH 
equation ( |3.13| ) possesses the following zero curvature representation 



Vx = 

Vt- 2 

which is equivalent to 



\-\m\ 

~~ 2 U ' JL 



-u - A" 1 
~mu\ + ju — tA _1 j 



(3.14) 



2 



y, 



ipt_ 2 = \u x ^) - wp x - A l ip x . 
Eq. (|3.15 ) is coinciding with the one in Ref. ]7). 



(3.15) 



In the cases of ai = 0, a 2 — 1, % = and ai = 0,a2 = 0,a3 = 1, i.e. GL] 



we can write them in a uniform expression: 

G-i = e ex , € = ±1. 



The first equation of Eq. fl3.5|) reads 

mt_i = (m, + 2em)e ea; (3.16) 

which is a linear PDE. 

Because G_ 2 = £ 1 ■ G-\ — [u + ei/ -1 )) e e:c , the second equation of Eq. (|3.5|) reads 

m t _ 2 = (m x (u + eM (_1) ) + 2m ( Ua . + 2eu + e ex (3.17) 

where m — u — u xx . This equation has the following zero curvature representation 



Vx = 

Vt-2 



1 

\-\m\ 
V- 2 y, 



y, 



(3.18) 
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where 

V-2 = -F(G_ 2 )A- 1 -y(G_i)A" 



2 



ex ( \{u x + 2eu + u^) + CX- 1 u + eut-V) + A" 1 

§^ + X7/.+ ^«/(- 1 ) + im(w + eu(- 1 ))A+iA- 1 -| (u x + 2eu + u^) + eA" 1 



Eq. (|3. 18 ) can be changed to the following Lax form 



ijj t _ 2 = (u + eu(-V + A" 1 ) e ex i/j x - \ (u x + 2eu + + A~ x ) e tx ifj. 1 ' ' 
Both of the two cases: e = ±1 for Eq. ( |3.17| ) are integrable. 

3.2 r-matrix structure for the Neumann- like CH system 

Consider the following matrix (called 'negative' Lax matrix) 



where 



A. (A) = -( P ,q)\-i + j2^, (3.21) 

j=i J 

N 2 

B- (A) = A^ + ^A- 1 -^^^-, (3.22) 

N ^2 



C-(A) = | A -2 - (p,p) A -1 + ^2 ~)~~~\~.- ( 3 - 23 ) 

j=l 3 

We calculate the determinant of £-(A): 



iA 2 detL_(A) = -iA 2 TrL 2 _(A) = -iA 2 (A 2 _(A) + i?_(A)C_(A)) 

1 N E~ 

= E^+Ex^- (^4) 

j=-2 3=1 
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where Tr stands for the trace of a matrix, and 

8 

H-i = \(p,p)-\(q,q), (3.25) 
#o = (p, q) (Ap, g) - (p, g) 2 , 
#i = - (p,g) (A _1 p,g) + (p,g) 2 , 

£7 = (p, g) Aj-pj-q,- - - ((g, g) Xj + \)y) 

1 

2 



\([p,p)\ 3 --^Jq 2 ] +\\]T-, 7 = 1,2,. ..,7V, (3.26) 



(PjQi-PiQj f 
Xj - Xi 



Let 



N 

F * = J2 Xk i +lE i> k=-l,-2,..., (3.27) 

3=1 



then it reads 



*i = i<A fe+ Vp)-i<A fe+1 g,g> 



i - 2 

+ 2 E « Ai+2 9> < Afc_i P.P> " < A '" +2 P> ^> < Afc_i P. 9» > ( 3 - 28 ) 



j=k 

k = -1,-2,-3,.... 

Obviously, F_i — H_i. 

Now, we consider the following symplectic submanifold in M? N 



M= (g,p)e 



t,2N 



F = l - ((Ag, g) - 1) = 0, G = (Ag,p> = > (3.29) 



and introduce the Dirac bracket on M 

{/, 9}d = {/, g} + ({/, g} - {/, G}{F, g}) (3.30) 

which is easily proven to be bilinear, skew-symmetric and satisfy the Jacobi identity. 
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In particular, the Hamiltonian system (H^i)z>: q x = {q,H-i}z), p x = {p,H-i}d 
on M reads 

[H- X ) D : { ll = b ~ 'i^AQ, (3.31) 
(Ag,p) = 0, (Ag,g) = l. 

We call this a Neumann-like system on M. Let 

m = l±ii^>, (3.32) 
2(A 2 g,g) 

yi = qj, V2 = p, A = Aj, j = l,...,N. (3.33) 
Then, the Neumann-like flow (H_i)d on M exactly becomes 

y x = U(m,X)y, y= (y 1 ,y 2 ) T , (3.34) 



which is nothing else but the CH spectral problem (|2.2 ) with the potential function m. 
Therefore, we can call the canonical Hamiltomnian system (|3.31|) the Neumann-like 
CH system on M. 

A long but direct computation leads to the following key equalities: 

M-(A), A-(jj)}o = {B_(A),B_(f,)} D = 0. 

(A 2 g,g) \/I A 

IA// (C_(AM_(/,)-C_(/iM_(A)), 



{A*q, q) 
2A/x 



{A-(\),B.(n)} D = -J-(B-.Qa) - fl_(A)) + y£_(Ai) + -S_(A) 

B_(A)B_(/x), 



(A 2 g, g) 

{A_(A),C_(/i)}D = ^(C_(A)-C_(/i))-|c_(/.)-V(A) 

(A 2 g, g) 2 (A 2 g, g) // 

{5_(A),C_0u)} o = -J—(A4fi)-A„(\)) + ^A_(fi) + -A4\) 

/i — A A yU 

IA// B (A).l (,,). 



(A 2 g, g) 
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Let L{ (A) = L_ (A)®/ 2 x2, L% (fi) = J 2 x2®£- (fx), where L_ (A) , L_ (li) are given 
through Eq. (|3.20|) . In the following, we search for a general 4x4 r- matrix structure 
rj~ 2 (A,/i) such that the fundamental Dirac-Poisson bracket: 



{L_ (A) f L_ (li)} d = [r u (A, Li) , (A)] - [r 2 l (/x, A) , L 2 (/x)] 
holds, where the entries of the 4x4 matrix {L_ (A) f L_ are 

{L_ (A) f L 0*)}^ = {L- (A) fcm , L_ , kj,m,n = 1, 2, 

and (A, ytx) = Pr± 2 {K A 4 ) ^ with 



- I ^2x2 



i=i 



a,- (X) cr,- 



/ 1 \ 

10 

10 

\ 1 ) 



(3.35) 



where a'j are the Pauli matrices. 



Theorem 3 



r i 2 0,/i) = 



2A 



/x(/x - A) 

zs an r-matrix structure satisfying Eq. j \3. Sty ), where 



p + s- 



s~ 



A(l + A(Ap,p)) (00 
1 



2/x (A 2 g, q) 
2A/x 



(A 2 q,q) 

( o 









1 







. A(l+4<Ap,p» n 



<8> 





1 



-B_(A) 

2A_(/x) B_(A) 

\ 





(3.36) 



Apparently, our r-matrix structure ( |3.36p is of 4 x 4 and is different from the one 
in Ref . |0| . 
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3.3 Integrability 

Because there is an r-matrix structure satisfying Eq. ( |3.35| ) 



{L 2 _(A)«L 2 (n)} D = [fr 2 (A )At ),Lr(A)] - [f 2 - x A) , (//)] (3.37) 

where 



1 i 

T7 ' ' 



(\ p) = e e^") 1- * w ) 1-1 m r « ^ ^ w ^ i m - 

A:=0 Z=0 

= 12, 21. 

Thus, 

4 {TrL 2 _ (A) , TiL 2 _ ( Ji )} D = Tr {L 2 _ (A) f L 2 _ (/x)} D = 0. (3.38) 
So, by Eq. Q3.24Q we immediately obtain the following theorem. 

Theorem 4 The following equalities 

{E~, EJ} D = 0, {H u EJ} D = 0, {F k , Ej} D = 0, (3.39) 
i,j = l,2,...,N, l = -2, -1,0,1, A: = -1,-2,..., 

hold. Hence, the Hamiltonian systems (Hi) D and {Fk) D on M 

(Ht) D : q x = {q,m} D , p x = {p,Hi} D , 1 = -2, -1,0,1, (3.40) 
(F k ) D : q tk = {q,F k } D , Ptk = {p,F k } D , k = -1, -2, . . . , (3.41) 

are completely integrable. 

In particular, we obtain the following results. 



Corollary 1 The Hamiltonian system (H-i)d defined by Eq. $3. 3\ ) is completely 
integrable. 

Corollary 2 All composition functions f (Hi, F k ), f E C°° (til), k = —1, —2, 
are completely integrable Hamiltonians on M. 



16 



3.4 Parametric solution of the negative order CH hierarchy 
restricted onto M 



In the following, we consider the relation between constraint and nonlinear equations 
in the negative order CH hierarchy (|3.5|). Let us start from the following setting 



N 

G\ = ^VA„ (3.42) 
i=i 

where G\. x = 1, and VAj = AjgJ is the functional gradient of the CH spectral problem 
( |2.2|) corresponding to the spectral parameter Xj (j — 1, N). 
Apparently Eq. ( |3.42j ) reads 

(Aq,q) = l. (3.43) 
After we do one time derivative with respect to x, we get 

(Ap,g> = 0, p=q x . (3.44) 



This equality together with Eq. ( |3.43| ) forms the symplectic submanifold M we need. 



Apparently, derivating Eq. ( |3.44|) leads to the constraint relation ( |3.32| ). 



Remark 1 Because M defined by Eq. ( \3. 2% is not the usual tangent bundle, i.e. 
_t TS n-i = ^ q p j e R 2V| F=(q,q)-1 = 0,G= (q,p) = oj and Eq. fl3J2]) is 



not the usual Neumann constraint on TS N 1 , Eq. ( 3. 33\ ) is therefore a new kind of 



Neumann system. In subsection 3.3 we have proven its integrability. 

Since the Hamiltonian flows (H_i)d and (F^d on M are completely integrable and 
their Poisson brackets {H-i,Fk}D — (k = —1,-2,...), their phase flows g x H _ v g l p k 
commute 0. Thus, we can define their compatible solution as follows: 

q ^ tk) )-a* ( q{x ° A) ) k--l-2 (3 45) 

where x°, t° k are the initial values of phase flows g x H _ x , gp . 
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Theorem 5 Let q(x,t k ), p(x,tk) be a solution of the compatible Hamiltonian systems 
(H-i)d and (Fk)n on M. Then 

= l + 4(Ap(x,t k ),p(x,t k )) 
m 2(A 2 q(x,t k ),q(x,t k )) 1 ' ' 

satisfies the constrained negative order CH hierarchy (i.e. constrained to M.) 

m tk = JC k+1 -l, k — — 1, —2, (3.47) 

where the operators J and C^ 1 are given by Eqs. I\2.b\ ) and (\S.b\) , respectively. 

Proof: On one hand, the recursion operator £ acts on Eq. ( |3.42j) and results in the 
following 

JC k+1 -G\ = J-(A k+2 q,q) 

= m x (A k+2 q, q) + Am (A k+2 q, p) 

= 2(1 + 4(A f 2 ' P)) ((A\q) (A k +\p) - (A\p) (A k +\q)) . 
\A 2 q, q) 

(3.48) 

In this procedure, Eqs. Q2.4|) and Q3.31|) are used. 

On the other hand, we directly make the derivative of Eq. ( |3.46|) with respect to 
t k . Then we obtain 

4 {A 2 q, q) (A 2 p, p tk ) - (1 + 4 {Ap, p)) (A 2 q, q tk ) 
(A 2 q,q) 2 

where q = q(x, t k ), p — p(x, t k ). But, 

Qt k = {q,F k } D , Pt k = {p,F k } D , k = -1, -2, (3.50) 
where F k are defined by Eq. (|3.28|), i.e. 



-i 

]T {(A J+2 q,q) A k ->p- (A> +2 q,p) A^q) , (3.51) 

j=k 

-1 

+ ^ « A J+2 9, p) A*-fy - (A J+2 p, p) A fe "^) . (3.52) 



Therefore after substituting them into Eq. ( |3.49|) and calculating it, we have 



m tk = 2(1 + 4<A g P)) ((A 2 g,g) (A k+ \P) ~ (A\P) (A k+ \q)) 
{A 2 q, q) 

which completes the proof. 

Lemma 1 Let q, p satisfy the integrable Hamiltonian system (H_i) D . Then on the 
symplectic submanifold M, we have 



(q,q)-4(p,p) = 0. (3.53) 

2. 

u= (q(x,tk),q(x,t k )) , k,= -1, -2, .., (3.54) 



satisfies the equation m = u—u xx , where m is given by Eq. ( $-4(\ ), and q(x, p(x 



be a solution of the compatible integrable Hamiltonian systems (H_i) D and (F k ) D 
on M. 



Proof: 



{(q,q) -4(p,p))x = 2(q,p) -8{p,p x ) 
= 



completes the proof of the first equality. 
As for the second one, we have 



u - u xx = (q, q) - 2 (p, p) - 2 (q, ^q - Knkq 

= g te>i) ~ 2 (p-'P) + m 

= m. 



In particular, with k = —2, we obtain the following corollary. 



19 



Corollary 3 Let g(x,t_ 2 ), p(x,t^ 2 ) be a solution of the compatible integrable Hamil 
tonian systems (H_i)d and (-F- 2 )d on M. Then 

l+4(Ap(x,t_ 2 ),p(x,t^ 2 )} 



m 
u 



mix, t- 



(3.55) 
(3.56) 



2(A 2 q(x,t_ 2 ),q(x,t_ 2 )) ' 
u(x, £_ 2 ) = (q{x, t- 2 ), q{x, £- 2 )) , 

satisfy the CH equation ( \3.1Q ) on M. Therefore, u = (q(x, t- 2 ), q(x, t- 2 )) , is a solution 
of the CH equation on M. Here U_i and F„ 2 are given by 

H-i = ^(p,p) -~(q,q) , 

ill 

F-2 = -(A~ 1 p,p)--(A~ 1 q,q) + -((q,q)(p,p)-(q,p) 2 ). 



Proof: Via some direct calculations, we obtain 

2(l + 4(Ap,p)) 2 w , 2 w , 

m ^ = 7-^ — — (\ A q>v w.p) - \ A ?>P) 

(A 2 g,g) 

And lemma 1 gives 

-J-u = -JC- l -G\ 

= J-<q,q> 

2(l + 4(Ap,p)) 



(A2g,g) 2 



«A 2 ^,g) (q,p) - (A 2 q,p) (q,q)) 



= m. 2 

which completes the proof. 



By Theorem |5|, the constraint relation given by Eq. Q3.32j ) is actually a solution of 
the constrained negative order CH hierarchy (|3.5| ) on M. Thus, we also call the sys- 
tem (H-i)d (i.e. Eq. ( |3.31D ) a negative order restricted CH flow (Neumann- 
like) of the spectral problem (|2.2| ) on the symplectic submanifold M. All 
Hamiltonian systems (F k ) D , k < 0,k 6 Z are therefore called the negative 
order integrable restricted flows (Neumann-type) on M. 

Remark 2 Of course, we can also consider the integrable Bargmann-like CH system 
associated with the positive order CH hierarchy Please see this in section 4- A 

systematic approach to generate new integrable negative order hierarchies can be seen 
in Ref. fiTT 
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3.5 Comparing parametric solution with peakons 

Let us now compare the Neumann-like CH system (H-i)d with the peakons dynamical 
system. 

Let Pj,Qj (j = 1,2, ...,N) be peakons dynamical variables of the CH equation 
( ^131) , then with § 



,t) = Y^Pj(t)e-\ x - Q '®\ (3.57) 



AT 

//(./'. 

i=i 

N 

m(x,t) = ^PfoWx-Qjit)), (3.5* 



where t = £_2 and 8{x) is the 5-function, the CH equation ( |3.13| ) yields a canonical 
peaked Hamiltonian system 



QM-M = ^ Pk ^ e ~ mt) ' Qm 

PAt) = ~M] = YLx Wag* {Q k {t) - QM e -W*(*)-«iWI 



(3.59) 

r iW = _ aoJ = ~ r o l^h=\ J-kWayn\wW) ~ WjK L )) c ^ J "~" 

with 

1 N 

H{t) = -^Pi(t)Pj(t)e-\ Qt ®- Q '®\. (3.60) 

In Eq. ( |3.59|) , "sgn" means the signal function. The same meaning happens in the 
following. 

A natual question is: what is the relationship between the peaked Hamiltonian 
system ( |3.59[ ) and the Neumann-like systems ( |3.31| ) and (F_ 2 )l Apparently, the peaked 
system (|3.59|) does not include the systems ( |3.31D and (F_ 2 ) because the system (|3.59|) 
is only concerned about the time part. 

By Corollary [| we know that 

N 

u = u(x,t) = (q(x,t),q(x,t)) = (3.61) 

3=1 

is a solution of the CH equation ( |3.13|) on M, where we set t_ 2 = t, and qj(x, t), pj = 
satisfies the two integrable commuted systems (H_i) D , (F_ 2 ) d on the symplectic 
submanifold M. 
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Assume Pj(t), Qj(t) are the solutions of the peaked system (|3.59|) , we make the 
following transformation (when Pj(t) < 0, we use a/ Pj (t) = i^J—Pj{t)) 



qj (x,t) = sJp~fc)e~^ x ~ Qm , 



which implies 



Pj(x,t) = -^qj{x,t) 



-sgn{x - Qj(t))qj. 



Hence, we obtain 



d d 2 1 

■^Pj(x,t) = -^qj{x,t) = -8(x - Qj(t))qj + -qj. 



(3.62) 



(3.63) 



(3.64) 



However, on M we have the constraints < Aq, q >= 1, < Aq,p >= which implies 

- 1 

Xjq^(x - Qj) = -, VxGl. (3.65) 

3=1 

This equality is obviously not true! So, the CH peakon system ( |3.59|) does not coincide 
with the nonlinearized CH spectral problem fl3.31| ). 

Let us now furthermore compute the derivative with respect to t. Inserting the 
peakon system ( ft.59| ), we get 



N 



nlj Pk ^ S9n ( X ~ ~ S 9 n (Qk(t) ~ Qj(t)) 



-\Q k {t)-Qj{t)\ 



. (3.66) 



k=l 



On the other hand, from the Neumann-like system {F-2)d we have 
qj(x,t) = {q j} F^ 2 } D 

= ^fpj - (q,p)vj 
i N 

= 2 qi ( ^ P ^) s 9< x - Qk(t))e^ x - Qk ^ - Xfsgnix - Qj{t))) . (3.67) 



fc=i 



Apparently, ( |3.66|) = (|3.67|) iff when x = Qj(x,t), i.e. for other x, they do not equal. 
Thus, the peakon system ( |3.59| ) is not the Neumann-like Hamiltonian system (F-2)d, 
either. 

So, by the above analysis, we conclude: the two solutions ( |3.57|) and (|3.61| ) 
of the CH equation ( |3.13| ) are not gauge equivalent. In the next subsection we 
will concretely solve Eq. ( [3.6 1|) on M in the form of Riemann-Theta functions. 
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3.6 Algebro-geometric solution of the CH equation on M 

Now, let us re-consider the Hamiltonian system [H-i)d on M under the substitution 
of A — > A _1 ,Aj — > Xj 1 (here we choose non-zero A, Xj). Then, the Lax matrix (|3.20|) 
has the following simple form: 



where 



L_(A) = -A 3 L CH (A), (3.68) 



Ach(X) Bcu(X) 
Cch(A) -A C h(A) 



(3.69) 



and the symplectic submanifold M becomes 
M = \(q,p) E R 2N 



F = 1((A- I g,g>-1)=0, G = (A^p) = 



(3.70) 



where A 1 = diag(X 1 , X N 



A direct calculation yields to the following theorem. 



Theorem 6 On the symplectic submanifold M the Hamiltonian system (H-i) d defined 



by Eq. ( ]3. 3f{) has the Lax representation: 



^L C h = [M C h,L ch ], (3.71) 



where 



4 4{A 2 q,q) 



M CK — ( 1 _ l+4(Ap,p) ^-i q ) • (3.72) 



Notice: On M the Hamiltonian system (H_i)d is 2A^ — 2-dimensional, that is, only 
there exist 2N—2 independent dynamical variables in all 2N variables q±, q^', Pi, ■■■Pn- 
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Without loss of generality, we assume the Hamiltonian system (H-\)d has the inde- 
pendent dynamical variables q%, ...,qN-i, Pi, ...,Pn-i {N > 1) on M. Then, on M we 
have 

N - 1 X 



q 2 N = - ~v: q2 v ( 3 - 73 ) 

3- 

N-l 



^ A? 

=1 3 

A; 



where gjv in the latter is given by the former in terms of q\, q^-i- 

We will concretely give the expression u = (q(x,t),q(x,t)) , t = t_ 2 in an explicit 
form. By Eq. (|3.73|) , We rewrite the entry Bqu (A) in the Lax matrix (|3.69|) as 

Sch(A) " "A^l 1 + ^ A - Xj q - 



3 



Q (A) 

K(xy 



(3.75) 



where 



AT-1 JV-l . . JV-l 

Q(X) = ^(A-A^ + ^^-^gJ II (A - A fe ) , (3.76) 

j=l j=l 3 k=l,k^j 

N 

K (A) = J](A-A Q ). (3.77) 

a=l 

Apparently, Q(X) is a iV — 1 (N > 1) order polynomial of A. Choosing its iV — 1 
distinct zero points fa, ... , fav-i, we have 

N-l 

q(a) = n( A -^)' ( 3 - 78 ) 

3=1 

N N-l 

(q,q) = E A «-E^'- ( 3 ' 79 ) 

a=l j=l 

Additionally, choosing A = Xj in Eqs. Q3.76Q and (|3.78|) leads to an explicit form of q 3 - 
in terms of fa: 

iV_1 A- 
q) = UjN JJ {Xj - fa), a jN = N 3 — , (3.80) 

1=1 llfe=ll A J - A k) 
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which is similar to the result in Ref. [0. By Eq. (|3.79[) , we get an identity about hi: 

N N-l N N-l 

Y a i N II ( A J ~ = Y Xa ~ Y 

3=1 1=1 a=l j=l 

Remark 3 The dynamical variable pj corresponding to q 3 - is 



N-l , N-l 



P.J 



dqj otjN ^A* 



dx 2<7, f--' cfc 

fc=i i=i,i^k 



therefore, 



N-l , iV-1 



9 J = a TfN-1 



* = r(E ^ II - *)) • ( 3 - 82 ) 



4ll/=i - fc =i — 1=1,1^* 

Substituting Eqs. ([3. 7Q) and (\3.8$\) into the Hamiltonian H_i = \{p,p) — §(<?,<?) 
directly gives an expression in terms of 

N N-l , N-l ,, , A - 1 , A 

V- 

(3.83) 



1 JV AT-1 , N-l „ 1 iV-1 1 AT 

E nf! fr^(E|f n + E*- E*. 

7=1 11/ = 1 l A J W fe=l l=ll^k 7 = 1 fe=l 



This is evidently different from the Hamiltonian function in Ref. (see t/iere section 3). 
Here our H_i comes from the nonlinearized CH spectral problemn, i.e. it is composing 
of a Neumann-like system on M, which is shown integrable in subsection 3.3 by r- 
matrix process. It is because of this difference that our parametric solution does not 
include the peakons (also see last subsection) and we will in the following procedure 
present a class of new algebro- geometric solution for the CH equation constrained on 
the symplectic submanifold M. 

Let 

7^ = AchM, j = l,..., N-l, (3.84) 
then it is easy to prove the following proposition. 

Proposition 1 

{m,Hj}D = {'Xi,'Kj}D = 0, {nj,Hi} D = Sij, i,j = 1,2, N- 1, (3.85) 
i. e. -Kj , fj,j are conjugated, and thus they are the variables which can be separated ^Zl 
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Write 



detL CH (A) = A 2 cll (A) + B Cil (A) C CH (A) 
1 (l " 



A 2 \ 4 z — f A - A c 

0=1 



(A a — Aat)A q j 



A(A-A iV )\4 ^ A-A c 
P(A) 



a^: (A) ' 



(3.86) 



where E a is defined by 



E a = -Pi+\ii-r a , r Q = x: ^-g^ 2 , (3 - 87) 

AT-1 N-l . . JV-1 

p(A) = ^n( A - A .)+E^x^ 2 n ( a - a ^)' ( 3 - 8g ) 

3=1 j=l - 1 k=l,ky£j 

and obviously P (A) is an N — 1 order polynomial of A whose first term's coefficient is 
7. Then we have 



1 



^=^-y j = l,...,N-l. (3.89) 
Notice: On M we always have ^2^=1 K^Ea = \- Therefore, we assume En,..., E^-i 



are independent. Then, 



N-l 

I 

E N 



4 ^ A a 

a=l 

Actually in Eq. ( |3.86| ) we used this fact. 
Now, we choose the generating function 

N-l 

W = J2 W i(^{Ea}Z=?) 

V - J m I p(A) 

.1=1 J K 
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where /xo is an arbitrarily given constant. Let us view E a (a = 1, . . . , N — 1) as actional 
variables, then angle-coordinates Q a are chosen as 

dW 

Q a = 777^, a = l,...,N -1, 



dE a 



i.e. 



^ ' 2v/AK(A)P(A) A Q (A — A Q )(A — Ajv) 



fc=l "'W 

\ \ N ~ 1 rv-k 

- ^V^E / (3-91) 

fc=l J W 

where 

oj q = -. dX, a — 1, . . . , N — 1. 

2v/AK(A)P(A) 

Therefore, on the symplectic submanifold (M 2N ~ 2 , dE a A dQ a ) the Hamiltonian 
function 

1 1 1 N 

0=1 

- (*<«> 

produces a linearized x-flow of the CH equation 



f) _ 9H-i _ Aq-Ajy 

^ a > x 8E a 2\ a ' 



E =0 

J -'a,x w ) 

as well as the Hamiltonian function 



(3.93) 



N 



F-2 = ^(Ap, P )-^(Aq,q) + ^{(q,q)(p,p)-(q,p) 2 ) =-^Ys X « E « 

a=l 

1 1 N-l > 2 , 2 

«=1 
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yields a linearized i-flow of the CH equation 

0F_ 2 



a,t 



Q 

Ea, t = 0. 



dE a 



\ 2 \ 2 
2X a ' 



The above two flows imply 

Qa 



Xn — A c 
2A„ 



x + (X N + X a )t - 2Q° a 



constant, a = 1, N — 1, 
where Q° a is an arbitrarily chosen constant. Therefore we have 

1 



N-l 

E 



(3.95) 

(3.96) 
(3.97) 

(3.98) 



UJa = -~ X + (Ajv + X a )t + Q°. 

fe=l J W> 

Choose a basic system of closed paths «j, (z = 1, . . . , iV — 1) of Riemann surface 
T: /j, 2 = XP (A) K (A) with N-l handles. Uj (j = 1, . . . , N — 1) are exactly N-l 
linearly independent holomorphic differentials of the first kind on the Riemann surface 
T. Let Qj be normalized as Uj = Yla=i r 3,1^1^ l - e - u j satisfy 



where B = (-Sy)(jv_i)x(iv-i) * s symmetric and the imaginary part ImB of B is a positive 
definite matrix. 



By Riemann Theorem [15| we know: fj,k satisfy 

N-l 



E 



fe=l 



JV-1 



0, (x, £) = rjJ Q? - - x + (Ajv + A,)t 



j = l,...,iV-l, 



iff //fc are the zero points of the Riemann-Theta function (P) = (A (P) — (ft — K) 
which has exactly iV — 1 zero points, where 



•A CP) 



wi, • • • , 



Pa 



P v 



= = (Mm) , • • ■ , <?W-i om)) t , 
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IK = (Ki, Kn^i) t G C^ -1 is the Riemann constant vector, Pq is an arbitrarily 
given point on Riemann surface T (O-function and the related properties are seen in 
Appendix). 



Because of 10 



±J\dlne(P) = C 1 (f), 



(3.99) 



where the constant C\ (F) has none thing to do with 0; 7 is the boundary of simple 
connected domain obtained through cutting the Riemann surface F along closed paths 
a>i, (3i. Thus, we have a key equality 

N-l 

(j, k = d (r) - Res A=001 Adln6 (P) - Res A=002 Adln 9 (P) , (3.100) 

k=l 

where ResA=oo fc (k = 1,2) mean the residue at points oo^: 

d 



(0, y/z-iP(z-l)K(z-l] 



2=0 



2=0 



00 2 = (0,- y/z^Piz-^Kiz- 1 ) 

Now, we need to calculate the above two residues. 
Lemma 2 

Res A=0Ol Adln9(P) = -2^- In 6 (0 + K + %) 

ox 

d 

ResA=oo 2 Adln6(P) = 2— In 6 (0 + K + r) 2 ) , 
where rji, r\2 are two different N — 1 dimensional constant vectors. 



(3.101) 
(3.102) 



Proof Consider the following smooth superelliptic curve F: // 2 = XP(X)K(X). 
Because AP (A) K (A) is a 2N-tYi order polynomials with respect to A, 00 is not a 
branch point, i.e. on F there are two infinity points 00 1 and oo 2 . All points P on 
f are denoted by (A, ±/i). On f we choose a group of basis of normalized closed 
paths a±, ■ ■ ■ , CKjv-i; Pi, • • • , Pn-i- They are mutually independent, and their intersec- 
tion number are 

at o ctj = Pi o Pj = 0, a;o Pj = 5ij. 

It is easy to see that Uj (j = 1, N — 1) are N — 1 linearly independent holomorphic 
differential forms on F. 
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Let us now come to calculate the residus of Xd In (P) at the two infinity points: 
ooi,oo 2 . At ooi, the j-th variable Ij of (z) produces the following result through 
multiplying by —1 (please note that the local coordinate at ooi and oo 2 is z = A -1 ): 



N-l „ z 

-Ij = fa + Kj + rju - r J,i / &i 

i=i Jo 



v 1 -n^U 

'3 "T- '\/ -r '/I,/ 



- /v ./ + >h-.i - £ '•/■/ / 

./0 



z 2 dz 



.o 2 V /AP(A)K(A) 

v-i 

+ ^ + '/I,, + E / ^ J= ; = , = — ^ — <i- 



7V-1 



~ Jo y/1 + 0(Z) 
N-l 

= 4> J + K J + r ] ^ + Y,r h iZ + 0{z 2 ) , 



where r} hj = JJ ± Uj, j = 1, • • • N - 1. 
Because 



5x 2 2-*> 2-*> dL rj ' 1 

3=1 1=1 3 



and (z) has the expansion formula 



N-l N-l ^ 



(z) = (0 + K + 770 - £ ^ W + O (^ 2 ) , 

j=l 1=1 3 



where 771 = (7/1,1, -, 77i,at_i) t . Therefore, 

<90 

(z) = ((f) + K + 771) - 2—z + O (z 2 ) . 
So, we obtain the following residue 

ResA=ooiAc?ln0 (P) = Res 2=O 2:~ 1 (jln0 (z) = Res z= Q- ^~ z — 

z (z) 
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1 -2G X + O (z) 
KeSz=0 z Q-2Q x z + O(z*) 
1 -2Q X + O (z) 
z -°zQ(l-2Q- 1 Q x z + 0(z 2 )) 

In a similar way, we can obtain the second residue formula. 



So, by Eq. ( |3.79[ ) and this Lemma, we immediately have 



(, <«, ( ) . t (z, t)> = £ *• - C. (f ) + 2f (in 9 + K + *> 



a=l 



&r V (0 + K + r)i) J ' 



(3.103) 



where the j-th component of rji (i — 1, 2) is r^j = J^ Wj. 

So, the CH equation ( 3.13|) on M has the following explicit solution, called the 
algebro-geometric solution 



u (x ,t) = R + 2l(^ + K + ^ 



dx\ e(0 + K + r/ 1 )y' 



(3.104) 



where R = J2a=i ~ @i (r) is a constant. 



Theorem 7 The algebro-geometric solution of the CH equation ( |3.13|) on M can be 
given through Eq. ( \3.104 )- 



Remark 4 Here the algebro-geometric solution ( \3. 1 04\ ) is smooth and occurs in the x- 
direction (spacial variable) derivative, and apparently differs from the piecewise smooth 
algebro-geometric solution in the t-direction derivatives given in Ref. J2|/. tt is also 
different from the results in Ref. [| ; [| ; becasue we are studying the CH equation 
constrained on M. In the paper, we do not need to calculate each qj (j = 1, N — 1) 



but the sum J2j=i ih which we know from Eqs. ( \3.7Q) and ( \3. 56] ). From the above 



subsection's comparison and these comments, therefore we think Eq. Q3.104 ) ^ s a class 
of new solution of the CH equation (\3.1ty) on M. Apparently, Eq. (\3. 1 OJ^j is more 
simpler in the form than in Ref. Ml. 
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4 Positive order CH hierarchy, integrable Bargmann 
system, and parametric solution 



4.1 Positive order CH hierarchy 

Let us now give the positive order hierarchy of Eq. (|2.2|) through employing the kernel 
element of the Lenard's operator J. 

Obviously, Go = cm~^ form all kernal elements of J, where c = c(t n ) G C°°(R) is 
an arbitrarily given fuctions with respect to the time variables t n (n > 0,n G Z), but 
independent of the spacial variable x. Gq and the recursion operator (|2.6|) yield the 



following hierarchy of the CH spectral problem (|2.2j ) 

m tk =cJC k -m-2, k = 0,1,2,..., (4.1) 

where the operators J and £ are defined by Eqs. ( |2.5|) and fl2lf ), respectively, and C 
and J have a further product form 

C = r l K = ^m-^d^m'^d- d 3 ) (4.2) 
J' 1 = im^d^m-s- (4.3) 



Because Eq. ( |4.1|) is related to the Camassa-Holm spectral problem ( |2.2|) for the case 
of k > 0, k G Z, it is called the positive order Camassa-Holm (CH) hierarchy. 
Eq. (|47T| ) has the following representative equations 

m to = 0, trivial case, (4.4) 
m tl = — c(m _ 2) +c(m~2) . (4.5) 



Apparently, with c = — 1 Eq. (fL5|) becomes a Dym type equation 

m tl = (m~2 j — (m~5 ) , (4-6) 



which has more an extra term — i^m I j than the usual Harry-Dym equation m tl 



_ i 
m 2 



. Therefore, Eq. ( |4. 1| ) gives an extended Dym hierarchy corresponding to 
the isospectral case: A t =0. 
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By Theorem 0, the whole positive order CH hierarchy (|4.1|) has the zero curvature 
representation 

U t „ - V Kx + [U, V k ] = 0, k > 0, k e Z (4.7) 

fc-i 

3=0 

where U is given by Eq. (|2.2|) , and Vj = V'(Crj) is given by Eq. ( |2.9| ) with G = Gj = 
D ■ Go = cO ■ m~^, j > 0, j E Z. Thus, all nonlinear equations in the positive CH 
hierarchy ( |4.1|) are integrable. Therefore we obtain the following theorem. 

Theorem 8 The positive order CH hierarchy ( U-A ) possesses the Lax pair 



1 

\-\m\ 



\y, 



_ Y-vjfe-i / — \Gj,x ~Gj \ \k-j ' 

2 vjr j?>;r 4 J 2 3 2 3fi 

k = 0,1,2,.... 



Eq. ( fO|) can be reduced to the following Lax pair 

IHh = Eto (1%^-^ - Gj\ k ~ j ip x ) , k — 0,1, 



(4.10) 



(4.11) 



In particular, the Dym type equation ( f4.6| ) has the Lax pair 

tpxx = \^~ \m\%p 

Remark 5 This Lax pair coincides with the one by the usual method of finite power 
expansion with respect to spectral parameter A. However, we here present the positive 
hierarchy ( \4-l\ ) mainly by the Lenard's operators pair satisfying Eq. Due to 

containing the spectral gradient VA in Eq. ( \2.4\ ), this procedure of generating evolution 
equations from a given spectral problem is called the spectral gradient method. 

In this method, how to determine a pair of Lenard's operators associated with the 
given spectral problem mainly depends on the concrete forms of spectral problems and 
spectral gradients, and some computational techniques. From this method, we have 
already derived the negative order CH hierarchy and the positive order CH hierarchy. 
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4.2 r-matrix structure and integrability for the Bargmann CH 
system 

Consider the following matrix (called 'positive' Lax matrix) 



where 



A+{x) = ^^m, (413) 

j=l 3 



j=l 



2(A<Z.9> ^ A - Aj 
We calculate the determinant of £+(A): 



idetL + (A) = \TrLl(X)^ 1 -(Al(X) + B + (X)C + (X)) 
N E + 



X-Xj 

3=1 J 



where Tr stands for the trace of a matrix, and 



1^,1=1 

Let 



^ Aj - A/ 



AT 

F fc = J]Aj£+, fc = 0,1,2,..., (4.19) 
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then it reads 
F k = - 



4(Ag,g> 2 



j=0 



(4.20) 



where fc = 0,l,2....A long but direct computation leads to the following key equalities: 

{A + (\),A + (f,)} = {F + (A),F + (/x)} = 0, 
{C + (A),C + (/x)} = (AA+(A) - M+M) , 



{A + (A),F + (/x)} 
{A+(\),C + (»)} 
{F + (A),C + (/x)} 



(Aq,q) 
2 



/i — A 
2 

/i — A 
4 

/x — A 



(jmB + (ji)-XB+(\)), 
(AC+(A) - aA.(ai)) 



(M+(m)-aa + (a)). 



(Ag, g) 



;AB+(A), 



Let (A) = L+ (A)<8) J 2X 2, ^2 (/-O = hx2®L + (/1), where L + (A) , L + (/x) are given 
through Eq. (|4.12|) . In the following, we search for a 4 x 4 r-matrix structure r± 2 (A, n) 



such that the fundamental Poisson bracket 13 



{L + (A) f L + (/x)} = [ r + (A, , L+ (A)] - [r+ (/x, A) , L+ (,,)] 
holds, where the entries of the 4x4 matrix {L + (A) ® L + (/i)} are 

{L + (A) f L + (/i)} W)mn = {L + (A) fem , L + (/xU , fc, I, m, n = 1, 2, 
and r 2 i (A, /x) = Pr i2 (A, /x) F, with 



'2x2 



i=i 



°j <» a i 



/ 1 \ 

10 

10 

\ 1 / 



(4.21) 



where are the Pauli matrices. 
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Theorem 9 

2A A 

(A, fi) = rrP + — — 2 S+ (4.22) 

fi{fi - A) (Ag, q) 

is an r-matrix structure satisfying Eq. ( \4-21\ ), where 

, \ / 

s + =*-®*-=i ! o r i o 



In fact, the r-matrix satisfying Eq. (|4.21|) is not unique |L9|. Obviously, this r- 
matrix structure is also of 4 x 4 and different from the one in Ref. 
Because there is an r-matrix structure satisfying Eq. (|4.21| ), 



{L\ (A) f L\ (/*)} = [f+ (A, „) , L+ (A)] - [f+ (/x, A) , L+ (/.)] (4.23) 



where 



i i 



V . ■ 



(A, A») = E E(^) 1 ^ ( A ) (A*) ^ ( A > /*) (^) fc ( A ) M > 

y = 12,21. 



fc=0 Z=0 



Thus, 



4 {Tr4 (A) , Tr^ (/,)} = Tr {4 (A) f L\ (/,)} = 0. 
So, by Eq. (|4.16|) we immediately obtain the following theorem. 

Theorem 10 The following equalities 

{E+,E+} = 0, {F k ,E+} = 0, 
i,j = l,2,...,N, fc = 0,l,2,..., 



(4.24) 



(4.25) 



hold. Hence, all Hamiltonian systems (F k ) 

QF k 

{F k ) : qt k = {q, F k } = -—, p tk = {p, F k } 



OF, 



dp ' ■ rtfe " J dq 

k = 0,1,2,... 



(4.26) 



are completely integrable. 
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Furthermore, we find the following Hamiltonian function 



is involutive with Ej~, Fk, i.e. 



{H + ,E+} = 0, {H + ,F k } = 0, (4.28) 
j = l,2,...,N, fc = 0,l,2,.... 



Here, £7+ are iV independent functions. 

Therefore, we obtain the following results. 

Corollary 4 The canonical Hamiltonian system (H + ): 

(* + >: (429) 

zs completely integrable. 

Corollary 5 yMi composition functions f(H + ,F^), f E C°°(WL), k = 0,1,2,..., are 
completely integrable Hamiltonians. 



Let 



m = — -j, (4.30) 
^ = Qj , X = \ j , j = l,...,N. (4.31) 



Then, the integrable flow (H + ) defined by Eq. ( |4.29| ) also exactly becomes the CH 



spectral problem (|2.1|) with the potential function m. 



Remark 6 Eq. ( \j.3C\ ) is a Bargmann constraint in the whole space M. . Therefore, 



the Hamiltonian system (H + ) is of integrable Bargmann type. 
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4.3 Parametric solution of the positive order CH hierarchy 

In the following, we shall consider the relation between constraint and nonlinear equa- 
tions in the positive order CH hierarchy (|4.1| ). Let us start from the following setting 

N 

G = -X>A„ (4.32) 
i=i 

where Go = —m~^ } and VAj = Aj-gJ is the functional gradient of the CH spectral 
problem (|2.2| ) corresponding to the spectral parameter Xj (j = 1, ...,N). 
Apparently Eq. ( |4.32j ) reads 

m = —^-2 (4.33) 
{Aq,q} 

which coincides with the constraint relation (|4.30| ). 

Since the Hamiltonian flows {H + ) and (F k ) are completely integrable and their 
Poisson brackets {H + ,F k } — (k — 0,1,2,...), their phase flows g^ + , g£ commute 
f|. Thus, we can define their compatible solution as follows: 



p(x,t k ) ) 9h+9f * \ p{x\t k 



9h + 9%\ IZo 2 I' *=0,1,2,... > (4.34) 



where x°, t° k are the initial values of phase flows g%+, g l p k - 

Theorem 11 Let q(x,t k ), p(x,t k ) be a solution of the compatible Hamiltonian systems 
(H + ) and (F k ). Then 

m = 2 (4.35) 

(Aq(x,t k ),q(x,t k )) 

satisfies the positive order CH equation 

m tk = -JL k -m--\ k = 0,1,2, (4.36) 



where the operators J and C are given by Eqs. ( \2.3j ) and ( \2.6j ), respectively. 
Proof: This proof is similar to the negative case. 

In particular, with k = 1, we obtain the following corollary. 
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Corollary 6 Let q(x,ti), p(x,ti) be a solution of the compatible integrable Hamilto- 
nian systems (H + ) and (i*i). Then 

m = m(x,ti) = 2, (4.37) 

(Ag(x,ti),g(x,ti)) 

is a solution of the Dym type equation ( U-(\ )- Here H + and F% are given by 

11 1 

H + = - (p,p) - - (q,q) 



2 \r>n g wn/ 4(Ag,g)' 

*i = -|^|| + ^((Ag,p) 2 -(Ag,g)(Ap,p)). 

By Theorem |Tl, the Bargmann constraint given by Eq. (|4.3Q| ) is actually a solution 
of the positive order CH hierarchy Q4.1|) . Thus, we also call the system (H + ) (i.e. 
Eq. ( [4.29| )) a positive order constrained CH flow (i.e. Bargmann type) of 
the spectral problem ( |2.2| ). All Hamiltonian systems (F k ), k > 0, k £ Z are 
therefore called the positive order integrable Bargmann flows in the whole 
R 2N . 

In a further procedure, we can get the algebro-geometric solutions for the positive 
order CH hierarchy. 
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Appendix 



Abel mapping and the ©-function 

1 If the genus of a Riemann surface is g, this surface is homomorphic to a sphere with g 
handles. Such a basic system of closed paths (or contours) a±, ...,a g ,f5\, ...,(3 g can be chosen 
such that the only intersections among them are those of ccj and $ with the same number i. 

Let the Riemann surface be covered with charts (Ui,Zi), where z% are local parameters in 
open domains Ui, the transition from Zi to Zj in intersections U fl Uj being holomorphic. If 
in any Ui a differential tpi (zj) dz{ with meromorphic (pi (zi) is given and in the common parts 
Ui fl f/j, (zj) dzi = ipj (zj) dzj, then we say that there is an Abel differential f2 on the whole 
surface with restrictions Q\ Vi = ipi (zi) dzi. The Abel differential is of the first kind if all the 
<Pi (zi) are holomorphic. There are exactly g linearly independent differentials of the first 
kind <jj\,...,u g . They are normed if j a ujj = 5^, which condition determines them uniquely. 
We shall always assume them normed. The numbers f j3 _ ujj = B^ are called /3-periods. The 
matrix B = (B^) has the following properties: 1) By = Bji, 2) r= Im B is a positive definite 
matrix. 

We consider a (/-dimensional vector A(P) = | j | , where Pq is a fixed point of the 
Riemann surface and P is an arbitrary point. This vector is not uniquely determined, but 
depends on the path of integration. If the latter is changed then a linear combination of a and 
/3-periods with integer coefficients can be added: (^4 (P))j ^ (A {P))j + Yli n i$ij + Yli m iBij, 
i.e. A(P) i-> A(P) + J2 n i$i + J2 m iBi, where 5i is the vector with coordinates Sij, Bi is 
the vector with coordinates Bij. Thus A(P) determines a mapping of the Riemann surface 
on the torus JJ = C s /T, where T is the lattice generated by 2g vectors {5i,B{\ (which are 
linearly independent over R). This mapping is called the Abel mapping, and the torus J is 
the Jacobi manifold or the Jacobian or the Riemann surface. 

The Abel mapping extends by linearity to the divisors: 

a (J2 n k p k) = Yl nkA • 

2 Abel Theorem. Those and only those divisors go to zero of the Jacobian by the Abel 
mapping which are principal. The latter means that they are divisors of zeros and poles of 
meromorphic functions on the surface. If P^ is a zero of the function, then > and is 
the degree of this zero. If Pk is a pole, then < and \nf~\ is the degree of this pole. 

Of special interest is the case of divisors of degree g with all n& = 1, i.e. of non-ordered sets 
of g points P\, P g of the Riemann surface. All the sets of such kind form the symmetrical 
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g-th power of the Riemann surface. The Abel mapping has the form 




The symmetrical g-th power is a complex manifold of the complex dimension g and it is 
mapped on the Jacobian which is a manifold of the same dimension. A problem of conversion 
arises (the Jacobi inverse problem). It can be solved with the help of the 6-function. 

3 For arbitrary P £ C 9 , let 

9 (P) = ex P i™ ( BZ > Z ) + 27ri ( P ' Z )} ' 




(BZ,Z) = ^2 BijZiZj, Z = (z 1 ,...,z g f 
a 

( P ,Z) = ^PiZii p = (Pi,-,P g ) T - 
i=i 

The series converges owing to the properties of the matrix B. The 0-function has the 
properties 

6(-P) = Q (P); 
0(P + 5 k ) = 0(P); 

Q{P + B k ) = Q(P)exp{-iri(B kk + 2 Pk )}. 
Note that the 6-function is not defined on the Jacobian because of the latter property. 

4 Riemann Theorem. There are constants K = = 1, ...,<? (Riemann constants) 

determined by the Riemann surface such that the set of points P\,...,P g is a solution of the 
system of equations 




= lj, L = {Ij} 6 J, j = 1, ...,g, 



if and only if P±, ...,P g are the zeros of the function B (P) = („4 (P) — L — K) (which has 
exactly g zeros). Note that while the function (P) is not uniquely determined on the 
Riemann surface (it is multivalued) its zeros are multivalued, since distinct branches of (P) 
differs by exponents. 
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5 We define now Abel differentials of the second and of the third kind. The Abel differential of 



z being the local parameter at this point. Such a differential is uniquely determined if it 
normed: f Viw = 0, Vi. 

The Abel differential of the third kind VLpq has only singularities which are simple poles 
at the points P and Q with the residues +1 and —1, respectively. It is uniquely determined 
by the same condition. 

6 Proposition. If z is a local parameter in a neighbourhood of the point P and u>i = fi (z) dz 
is the Abel differential of the first kind, then 



the second kind, Op , k = 1, 2, has the only singularity at the point P which is a pole of the 
order k+1. The differential can be represented at this point as dzT (holomorphic differential), 




1 d*- 1 
(k - 1)1 dz k ~ x 



z=0 



* = !,... 



and 




which is also seen in Ref. [11|. 
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